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Abstract 

A search for heavy leptons decaying to a Z boson and an electron or a muon is presented. 
The search is based on pp collision data taken at = 8 TeV by the ATLAS experiment at 
the CERN Large Hadron Collider, corresponding to an integrated luminosity of 20.3 fb“h 
Three high-transverse-momentum electrons or muons are selected, with two of them re¬ 
quired to be consistent with originating from a Z boson decay. No significant excess above 
Standard Model background predictions is observed, and 95% confidence level limifs on fhe 
producfion cross secfion of high-mass frilepfon resonances are derived. The resulfs are in- 
ferprefed in fhe confexf of vecfor-like lepfon and fype-III seesaw models. For fhe vecfor-like 
lepfon model, mosf heavy lepfon mass values in fhe range 114-176 GeV are excluded. For 
fhe fype-III seesaw model, mosf mass values in fhe range 100-468 GeV are excluded. 
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1 Introduction 

Searches for new particles often utilize decays to electrons or muons, which allow for full four-momentum 
reconstruction with a mass resolution better than that achievable using hadronic or semileptonic decay 
modes. Dilepton resonance searches led to the discovery of the 7/i/r [1, 2], the T [3], and the Z boson [4, 
5], and they have been used at the CERN Large Hadron Collider (LHC) to place strong constraints on a 
variety of new particles such as additional gauge bosons [6,7]. Searches for low-mass trilepton resonances 
have been used to constrain lepton flavour violation in muon and r lepton decays [8, 9]. 

High-mass trilepton resonances are motivated by several extensions of the Standard Model (SM). Vector¬ 
like leptons (VLL) are invoked to explain the mass hierarchy between the different lepton generations [10]. 
They also arise in composite Higgs models [11,12] and models of warped extra dimensions [13,14]. Such 
leptons have masses much larger than those of the SM leptons, and are defined as colourless, spin-1/2 fer¬ 
mions whose left- and right-handed chiral components have the same transformation properties under 
the weak-isospin SU(2) gauge group. Another set of models predicting trilepton resonances is based 
on the type-III seesaw mechanism [15], which explains the origin of small neutrino masses through the 
introduction of heavy SU(2) triplets with zero hypercharge. 

This article presents a search for high-mass trilepton resonances with the ATLAS detector, using a data 
sample corresponding to 20.3 fb“' of integrated luminosity collected in pp collisions at ^[s - ^ TeV at 
the LHC. This search uses data events with at least three charged leptons (electrons or muons), two of 
which are consistent with originating from a Z-boson decay. Several signal regions are defined to be 
sensitive to the pair-production of heavy leptons that decay to SM leptons and W, Z, or H bosons. The 
backgrounds, dominated by SM diboson production, are estimated using Monte Carlo (MC) simulation 
and control regions in data, and the predictions are validated in dedicated data samples. 

The results of the search are interpreted in the context of vector-like lepton [16] and type-III seesaw [17] 
scenarios in which the new heavy leptons decay through mixing with electrons or muons {{) induced by 
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off-diagonal Yukawa couplings. In the type-III seesaw model, the masses of the three heavy leptons are 
assumed to be identical. Feynman diagrams of the production and decay of the heavy leptons in both 
models are shown in figure 1. The heavy leptons are produced in pairs through Drell-Yan processes, 
with cross sections of roughly 34 fb and 844 fb for the VLL and type-III seesaw models, respectively, 
assuming heavy lepton masses of 200 GeV. The difference in the production cross section is due to the 
different gauge couplings of the models, as well as the additional neutral fermion in the type-III seesaw 
model. The heavy leptons decay via the mixing terms into an SM lepton and a W, Z, or H boson. The 
charged states L- exist in both models, and have decay modes to Wv, Zt, and H(-, the neutral state 
is only present in the type-III seesaw model, and decays to W£, Zv, and Hv. The charged lepton 
branching fractions approach S(L- ^ W-v) - 50%, S{L- Z£-) = 25% and S(L* ^ HI-) - 25% 
for mi± » mu, in accordance with the Goldstone boson equivalence theorem [18]; at lower masses, the 
branching fractions to H and Z decrease as they become kinematically disfavoured. For the neutral lepton, 
the branching fractions to 'Wl, Zv, and Hv are identical to those of the charged leptons to Wv, Zl, and 
H£, respectively. 



Figure 1: Feynman diagrams for the production and decay of new heavy leptons (L^, N^) to final states resulting 
in a trilepton resonance. Diagram (a) shows the pair production of two charged heavy leptons, and (b) shows the 
associated production of a charged and a neutral heavy lepton. 


Searches for heavy leptons were previously performed at LEP, excluding vector-like leptons with masses 
below mi± = 101.2 GeV using the L- W-v decay mode [19]. A search for type-III seesaw heavy 
leptons was performed by CMS in pp collision data at = 7 TeV, using non-resonant trilepton signa¬ 
tures to exclude seesaw fermions with masses below mL± = 180-200 GeV, depending on the branching 
fractions assumed [20]. 


2 The ATLAS detector 


The ATLAS detector [21] is a multi-purpose detector covering nearly the full solid angle ^ around the pp 
interaction region. The beam pipe is surrounded by the inner detector (ID), consisting of silicon pixel and 
microstrip detectors and a transition radiation tracker. The ID is enclosed in a superconducting solenoid 
providing a 2 T axial magnetic field, and performs charged particle tracking for \t]\ < 2.5. 


* ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector 
and the z-axis along the beam pipe. The z-axis points from the IP to the centre of the LHC ring, and the (/-axis points 
upward. Cylindrical coordinates (r, (p) are used in the transverse plane, (p being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 0 as ?/ = - In tan(0/2). 
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The calorimeter system surrounds the solenoid, and consists of electromagnetic and hadronic components. 
The electromagnetic calorimeter is a lead/liquid argon (LAr) sampling calorimeter, and comprises a barrel 
(|? 7 l < 1.475) and two endcaps (1.375 < I 77 I < 3.2). In the range \q\ < 2.5, the detector is finely segmented 
in 77 to provide good spatial resolution. The hadronic calorimeter (HCAL) uses steel/scintillator tiles in 
the barrel (It/I < 1.7) and copper/LAr in the endcaps (1.5 < \r}\ < 3.2). In the forward region (3.1 < 
It/I < 4.9), electromagnetic and hadronic calorimetry is performed using copper/LAr and tungsten/LAr 
technology. 

The muon spectrometer (MS) features high-precision tracking chambers interleaved with dedicated trig¬ 
ger chambers located in a toroidal magnetic field. The magnefic field is generafed by a sysfem of fhree 
large superconducfing air-core foroid magnefs, wifh a bending integral of abouf 2.5 T-m in fhe barrel 
and up fo 6 T-m in fhe endcaps. The precision fracking is provided by monitored driff fubes (|t/| < 2.7), 
complemenfed by cafhode strip chambers in the forward region (2 < |t/| < 2.7). Triggering is per¬ 
formed by resistive plate chambers in the barrel (|t/| < 1.05) and thin gap chambers in the endcaps 
(1.05 < |t/| < 2.4). 

Events are recorded using a three-level trigger system. The first level, implemented in hardware, reduces 
the event rate to less than 75 kHz using a subset of the detector information. The second and third 
levels are implemented in software, and reduce the event rate to less than 400 Hz using the full detector 
information. 


3 Object reconstruction and event selection 

The data were collected during 2012 using triggers requiring either an electron or a muon with transverse 
momentum relative to the beam axis, p-i, greater than 24 GeV. The triggered electron or muon must 
also satisfy loose isolation requirements. These triggers are supplemented by triggers without isolation 
requirements, but with higher px thresholds of 60 (36) GeV for electrons (muons). Only data taken while 
the ID, calorimeters, and MS were functioning normally are considered. Events are required to have a 
reconstructed primary vertex having at least three associated tracks with pj > 400 MeV, consistent with 
the beamspot envelope. If more than one such vertex is found, the vertex with the largest ^ of its 
associated tracks is chosen as the hard-scatter primary vertex. 

Electron candidates are selected as energy clusters within a small window of size At/ x A 0 = 0.075 x 
0.125 in the electromagnetic calorimeter matched to a track in the ID. They are required to fulfil fighf 
idenlificafion criteria [22], have |T/ciusterl < 2.47, and nol be in fhe fransifion region befween fhe barrel and 
fhe endcap calorimefer (1.37 < |T/ciusterl < 1-52), where T/duster is the pseudorapidity of the barycentre of 
the energy cluster. Muon candidates are selected as tracks reconstructed in the MS matched to tracks in 
the ID [23] and are required to satisfy |t/| < 2.5. The muon momentum is determined from combining the 
information from the two tracks. Muons and electrons are required to have transverse momenta greater 
than 15 GeV and to be isolated from tracks and calorimeter energy deposits using the criteria described in 
ref. [24]. To ensure that the lepton track is consistent with originating from the primary event vertex, the 
ID track is required to satisfy Ir/o/cr^Tol < ^ l^o sin 6\ < 0.5 mm, where do and zo are the transverse and 
longitudinal impact parameters of the track with respect to the primary vertex, respectively, and is the 
uncertainty on the transverse impact parameter. In order to ensure constant trigger efficiency as a function 
of lepfon px, af leasf one elecfron or muon musf have px > 26 GeV and a separation AR = ^|{Krfp■ + (A 0)2 
less fhan 0.2 from fhe friggered elecfron or muon. The frigger efficiency is evaluafed fo be larger fhan 
95% when all offline selecfion criferia are applied. 
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Jets are reconstructed from topological clusters built from energy deposits in calorimeter cells using the 
anti-^f jet algorithm [25] with a radius parameter of /? = 0.4. The measured jet energy is calibrated using 
Pi- and //-dependent corrections for instrumental effects (e.g. passive material and non-compensating 
response of the calorimeters) derived from MC simulations and in situ techniques applied to data, and is 
corrected for additional pp interactions per bunch crossing (pileup) [26, 27]. 

After energy calibration, jets are required to have pj > 30 GeV and \r]\ < 4.5. In order to suppress jets 
from pileup interactions, if a jet has pi < 50 GeV and \r]\ < 2.5, then at least 50% of the scalar sum of the 
Pi of all tracks associated with the jet is required to come from tracks associated with the primary event 
vertex. Jets are also required to pass jet-quality selections to reject jets reconstructed from non-collision 
signals, such as beam-related background, cosmic rays or detector noise. 

Since leptons and jet candidates can be reconstructed as multiple objects, the overlap between the various 
objects is resolved by applying the following procedure. If two electrons are separated by AT? < 0.1, 
the electron with the lower pi is removed. If an electron and a jet are separated by AT? < 0.2, the jet is 
removed. If an electron and a jet satisfy 0.2 < AT? < 0.4, and the jet’s transverse momentum also satisfies 
T’Tjet > 30 GeV -I- 0.05//T,e^ the electron is removed. If a muon and an electron satisfy AT? <0.1, the 
electron is removed. If a muon and a jet are separated by AT? <0.1, the jet is removed if its transverse 
momentum satisfies pijet < 0.5pT,/i if Pt,ij < 200 GeV, or PT,jet < 100 GeV if pi^^ > 200 GeV. Finally, 
if a muon and a jet not removed by the previous requirement are separated by AT? < 0.3, the muon is 
removed. 

The missing transverse momentum, and its magnitude, are calculated from the vector sum of 
the transverse momenta of all calibrated electrons, muons, r leptons, jets, and all topological calorimeter 
clusters of energy not associated with other objects with |? 7 | < 4.9. 

Events are required to have at least three leptons (electrons or muons) passing the selection requirements 
above. At least one pair of leptons with the same flavour and opposite electric charge must have an 
invariant mass within 10 GeV of the Z boson mass, mz [28]. Events with four leptons consistent with the 
decay of two Z bosons, also within 10 GeV of mz, are vetoed. Eor the remaining events with four leptons, 
the lepton closest in AT? to the Z boson candidate, referred to here as the off-Z lepton, is used to form 
the trilepton mass. Eor the range of heavy lepton masses considered in this analysis, the Z boson and the 
off-Z lepton tend to be collimated; hence, to improve the signal to background ratio, events where the Z 
candidate and the off-Z lepton are separated by AT? > 3 are vetoed. 

Eor simulated events with an Lr decaying to three leptons with pi > 15 GeV and \r]\ < 2.5, of which two 
originate from a Z boson and have an invariant mass within 10 GeV of mz, the efficiency of this event 
selection for the Z + e {Z + p) decay channel ranges from 20% (36%) at mi± - 100 GeV to 35% (38%) at 
mi^± - 400 GeV. The determination of the efficiency is discussed in section 8. 

Since the heavy leptons are produced in pairs, in addition to the identified L- Z + i decay, signal 
events contain either a second or an V°, which decays to a W, Z, or H boson and a charged or neutral 
lepton. A large fraction of events therefore contain a fourth lepton and/or a hadronically decaying boson. 
The sensitivity of the analysis is significantly improved by separating the events selected above into the 
following three exclusive categories: 

• at least four leptons are required using the same identification criteria as described above. 

• 'i{+]y. exactly three leptons are required, along with two jets with an invariant mass satisfying 
mw - 20 GeV < mjj <150 GeV, where m^ is the W boson mass [28]. 
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• 3^-only: the event does not fulfil the criteria of either the Ai or the Sf+jj categories. 

Subdividing the category based on the presence of a dijet does not significantly improve the sensitivity 
due to the small number of expected events with both a fourth lepton and a hadronically decaying boson. 
Finally, events are separated into two channels based on whether the off-Z lepton is an electron or a muon. 
This classification results in six independent signal regions. 

The search is performed by looking for a narrowly peaked excess of events in the distributions of the mass 
difference defined by Am = mj,e - where the invariant mass of the two leptons associated with the 
Z-boson decay is subtracted from the trilepton invariant mass. This reduces the impact of the lepton 
momentum resolution, and thus enhances the narrow resonance structure of the signal. The resulting 
reconstructed width in Am is 5.9 GeV (15.5 GeV) for a mass hypothesis of mL± = 120 GeV (400 GeV), 
while the corresponding width in the trilepton invariant mass is 7.3 GeV (18 GeV), for final states where 
the off-Z lepton is an electron. For final states where the off-Z lepton is a muon, the corresponding width in 
Am is 5.1 GeV (31.5 GeV) for a mass hypothesis of mi± - 120 GeV (400 GeV), while the corresponding 
width in the trilepton invariant mass is 6.7 GeV (33.5 GeV). The intrinsic width of the resonance is a few 
MeV at mi± = 120 GeV, rising to 0.5 GeV at mi± = 400 GeV. 


4 Monte Carlo simulation 

The analysis uses MC samples of VLL and type-III seesaw events generated with MadGraph 4.5.2 and 
5.2.2.1 [29], respectively, using the CTEQ6L1 [30] parton distribution functions (PDF) and the AU2 
underlying event tune [31]. Showering is performed with Pythia 8 [32]. Decays of the heavy leptons 
in the VLL model are performed using bridge [33], while decays in the type-III seesaw samples are 
performed by MadGraph. Lor the type-III seesaw model, the charged and neutral heavy leptons are 
generated with identical masses. Vector-like lepton samples are generated for eleven mass hypotheses for 
100 GeV < mi± < 400 GeV, while the type-III seesaw samples are generated for ten mass hypotheses for 
100 GeV < < 500 GeV. The cross sections for both samples are calculated at leading order (LO) 

in QCD. 

The main backgrounds originate from SM diboson production, in particular WZ and ZZ production. 
Contributions from WZ {ZZ) are modelled using the Sherpa [34] MC generator version 1.4.3 (1.4.5), 
using the internal showering algorithm [35-37], with the CTIO [38] PDL set and normalized to the next- 
to-leading-order (NLO) prediction from vbfnlo-2.6.2 [39]. The generation includes up to three additional 
parton emissions in the matrix element. Samples of simulated events based on the NLO generator powheg- 
Box [40] are used to derive systematic uncertainties on the shapes of distributions predicted by Sherpa. 
The diboson samples are showered with Pythia 8, and use the CTIO PDL set and AU2 underlying event 
tune. 

Drell-Yan production in association with a photon that converts in the detector, denoted Z -i-y, is modelled 
using Sherpa 1.4.1, also using the CTIO PDL set and including up to three additional parton emissions 
in the matrix element. Production of top-quark pairs in association with a IT or Z boson {tt + V) and 
triboson production (VVV^*^) are modelled using MadGraph 5.1.3.33, with Pythia 6.426 for the parton 
shower and hadronization, AUET2B underlying event tune [41], and the CTEQ6L1 PDL set. The tt + V 
processes are normalized to the corresponding NLO cross sections [42, 43], while the Z -i- y and wy(*) 
processes are normalized to their LO cross sections from the respective generator. 
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Table 1: Summary of the primary signal and background MC samples used in this analysis. The generator, parton 
shower and hadronization, PDF, underlying event tune, and the order of the cross-section calculation are shown for 
each sample. 


Process 

Generafor 

Parton shower and hadr. 

PDF sef 

UE fune 

Cross section 

VLL 

MadGraph 4.5.2 

Pythia 8 

CTEQ6L1 

AU2 

EG 

Seesaw 

MadGraph 5.2.2.1 

Pythia 8 

CTEQ6L1 

AU2 

EO 

WZ 

Sherpa1.4.3 

Sherpa 

CTIO 

Sherpa 

NEO 

zz 

Sherpa1.4.5 

Sherpa 

CTIO 

Sherpa 

NEO 

tt + Wiz 

MadGraph 5.1.3.33 

Pythia 6.426 

CTEQ6E1 

AUET2B 

NEO 

yyy(*) 

MadGraph 5.1.3.33 

Pythia 6.426 

CTEQ6E1 

AUET2B 

EO 

Z + y 

Sherpa 

Sherpa 

CTIO 

Sherpa 

EO 


For all samples, the response of the ATLAS detector is modelled using the geant4 toolkit [44, 45]. Addi¬ 
tional pp interactions in the same or nearby bunch crossings are included in the simulation by overlaying 
minimum-bias interactions modelled with Pythia 6.425 onto the hard-scatter event. The simulated events 
are reweighted to reproduce the distribution of the average number oi pp interactions per crossing ob¬ 
served in data. The generator, parton shower, PDF set, underlying event tune, and accuracy of theoretical 
cross section for the primary MC samples used are summarized in table 1 . 


5 Background estimation 


Standard Model processes containing three or more lepton candidates can be classified into two categor¬ 
ies. The first category consists of events with three prompt leptons produced in the decays of electroweak 
gauge bosons, which are estimated using the simulated samples described above. The second consists 
of events where at least one reconstructed lepton arises from a misidentified jet, hadron decay, or photon 
conversion, and is referred to as reducible background. For muons, reducible backgrounds arise from 
semileptonic b- or c-hadron decays and from in-flight decays of pions or kaons. Reducible electron 
backgrounds can arise from semi-leptonic b- or c-hadron decays, photon conversions and misidentified 
hadrons or jefs. Drell-Yan producfion of a lepfon pair wifh an associafed photon fhaf converts in fhe 
defecfor and is reconstrucfed as an isolafed lepfon (Z -i- y) is esfimafed using simulafion. The remainder of 
the reducible background is estimated by scaling control samples in data, following a method similar to 
that described in ref. [46]. The control samples consist of events with one or more leptons that do not sat¬ 
isfy the nominal selection criteria, but instead satisfy a set of relaxed criteria, defined separafely for each 
lepfon flavour. The evenfs are weighfed with scale factors computed for each such lepton, defined as fhe 
rafio of misidenfified or non-prompf lepfon candidafes fhaf satisfy fhe nominal criferia to those which only 
fulfil fhe relaxed criferia. For elecfrons, fhe idenfificafion requiremenf is changed from fighf fo loose [22]. 
For muons, fhe requiremenfs on fhe lepfon isolafion and on Ir/o/f^/ol are relaxed. The scale faclors are 
measured as a funcfion of fhe candidate’s pj and rj in samples of dafa fhaf are enriched in non-prompf 
and misidenfified lepfon candidafes. The contamination from prompt leptons in the background-enriched 
samples is removed using simulation. 

The background estimates are validated in four validation regions. The high-A/? region consists of events 
where the Z boson candidate and the off-Z lepton are separated by AT? > 3. The background composition 
in this region is similar to that in the signal regions. The off-Z region confains evenfs with exactly three 
leptons, where no opposite-sign same-flavour pair of leptons is reconsfrucfed wifh an invarianf mass 
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within 20 GeV of mz- This region is designed to test the Z+y background estimate. The ZZ region consists 
of events with two reconstructed Z boson candidates. These three regions are completely orthogonal to 
the signal selection. The fourth region is designed to test the modelling of SM WZ production. Events in 
this region must have exactly three leptons, zero jets, 40 < < 100 GeV, and 40 < < 90 GeV, 

where - cos(A(;i>)) is the transverse mass and Acp is the azimuthal angle between the 

missing transverse momentum and the off-Z lepton with momentum p^. This region is not completely 
orthogonal to the signal regions, but signal processes are expected to account for less than 3% of the 
expected event yield for type-III seesaw leptons with mL± >160 GeV. 

The expected and observed numbers of events are given in table 2 for all validation regions, separately 
for the Z + e and the Z + p channels. The largest difference is seen in the off-Z region in the Z + p flavour 
channel, where there is a deficit in the data corresponding to 2.3 standard deviations (cr). The region is 
dominated by contributions from ZZ, where only three leptons pass the selection requirements and no 
same-flavour, opposite-sign lepton pair is reconstructed with invariant mass within 20 GeV of mz- In the 
other seven regions, agreement better than 1.5cr is observed. Figure 2 shows the Am distributions for the 
high-A/? and ZZ validation regions in the Z + e and the Z + p flavour channels. 


Table 2; Summary of the number of events observed and predicted for each validation region. The uncertainty on 
the background prediction is the total systematic uncertainty. The difference between the observed and predicted 
number of events divided by the combined statistical and systematic uncertainty on the prediction is also shown. 


Channel 

Validation Region 

Data 

Background Prediction 

Data-Bkgd 

t^bkgd 

2 + 6 

High-AR 

239 

239 ± 14 

0.0 

2 + 6 

Off-Z 

360 

349 ±44 

±0.2 

2 + 6 

ZZ 

39 

37 ± 2 

±0.3 

2 + 6 

WZ 

140 

133 ± 10 

±0.4 

Z + p 

High-AR 

302 

301 ± 12 

±0.1 

Z + p 

Off-Z 

163 

200 ± 8 

-2.3 

Z + p 

ZZ 

74 

63 ± 3 

±1.2 

Z + p 

WZ 

222 

193 ± 14 

±1.5 


6 Systematic uncertainties 


Systematic uncertainties are assigned to the signal and background predictions derived from simulation 
to account for possible modelling inaccuracies. Sources of uncertainty affecting all simulated signal and 
background processes are the cross sections of SM processes, trigger efficiencies, lepton energy scales 
and resolutions (LES/LER), jet energy scale and resolution (JES/JER), lepton reconstruction and selec¬ 
tion efficiencies, MC statistical uncertainties, and luminosity. The cross-section uncertainties, evaluated 
for the SM background samples, include renormalization and factorization scale and PDF uncertainties. 
The scale uncertainties are determined by varying the renormalization and factorization scales up and 
down by factors of two. PDF uncertainties are obtained using the PDF4EHC working group recommend¬ 
ations [47]. Scale and PDF uncertainties are added in quadrature to obtain the total uncertainty on the 
inclusive cross section. For the dominant WZ and ZZ backgrounds, the resulting theoretical uncertainty 
on the NFO predictions from vbfnlo are 7.6% and 4.3%, respectively. A further uncertainty is assigned 
to the WZ and ZZ backgrounds to account for potential generator-level mismodelling of the shape of the 
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Figure 2: The Am = - m(*(- distributions for the high-A/? validation region (a) and (b) and the ZZ validation 

region (c) and (d). The figures (a) and (c) show the distributions for the Z + e final states, while (b) and (d) show the 
Z+jj final states. The error bars on the data points represent statistical uncertainties, and the shaded band represents 
the systematic uncertainty on the background prediction. 
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Am spectrum. The uncertainty is the difference between the predictions from Sherpa and powheg-box, 
symmetrized around the value from Sherpa. 

For the Z + y backgrounds, an uncertainty of 30% is assigned to the modelling of prompt photons con¬ 
verting in the inner detector, based on comparisons of conversion processes 'mZ ^ ee events between 
data and simulation. The reducible backgrounds are assigned an uncertainty related to the data-driven 
scaling procedure described in section 5, primarily due to the extrapolation of the scale factors from the 
measurement sample to the signal regions and to the correction for the presence of prompt leptons in 
the background-enriched samples. The uncertainties on the electron factors range from 24% to 30% as a 
function of pT> and the uncertainties on the muon factors range from 25% to 50%. 

The uncertainties on the lepton reconstruction and selection efficiencies, energy scales, and energy res¬ 
olutions [23, 48, 49] affect all simulated backgrounds, with combined uncertainties of 1% to 2% on the 
normalizations. The jet energy scale and resolution uncertainties [27, 50] only significantly affect the 
37’-Pjj signal regions, with a total uncertainty of 3%. Statistical uncertainties due to the finite number of 
events in the MC samples range from 1% to 5%. The luminosity uncertainty is 2.8%, and is derived 
using the same methodology as that described in ref. [51]. In total, the systematic uncertainties on the 
background predictions in each signal region range from 6% to 9%. 

The largest sources of uncertainty affecting the signal predictions are the lepton reconstruction and se¬ 
lection efficiencies, the luminosity, and, for the 3{+]] category, the jet energy scale and resolution. The 
total uncertainties on the signal normalizations range from 3% to 7% depending on the signal region and 
«IL±- 


7 Signal and background model 

The numbers of signal and background candidate events in data are determined from an unbinned maximum- 
likelihood fit of a combination of signal and background models to the Am distributions in each signal 
region. The details of the fit procedure and the models are described below. 

The signal and background processes are modelled with probability density functions (p.d.f.s). The para¬ 
meters of the p.d.f.s are determined from fits to the background estimates described in section 5. The 
fit to data using the combined signal and background model is performed simultaneously on the three 
categories for each of the two flavour channels. In each signal region, the normalization of the dominant 
background (ZZ or ITZ) is a free parameter in the fit. The normalizations of all other backgrounds are 
constrained to fluctuate according to Gaussian probability distributions with mean and width values equal 
to the estimates and the total uncertainties before fitting. The uncertainties on the p.d.f. parameters are 
incorporated as Gaussian-distributed nuisance parameters. 

The VLL and type-III seesaw signal models are parameterized separately as the sum of a Voigtian function 
(the convolution of a Breit-Wigner and a Gaussian function) for the trilepton resonance peak and a Landau 
distribution for the combinatorial part of the signal, where the three reconstructed leptons do not originate 
from the same Lr decay. The signal parameterization at a certain heavy lepton mass mL± is given as a 
function of Am by the following expression: 


S{mL±) = fyFy{Am-, Ty, my, cry) + (1 - fy)F]_{Am\ m^, ctl), (1) 
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where /v denotes the fraction of events in the resonance peak (Voigtian function), Fy, tn\ and cry are the 
width, mean, and Gaussian smearing terms of the Voigtian function Fy, and m\^ and ctl are the parameters 
of the Landau distribution Fl- The six parameters are determined at each simulated mass point by fitting 
S{mL±) to the simulated Am distributions, separately for the two flavours. For mass points mL± that lie 
between those assumed in the MC samples, the parameters of the signal templates are obtained by linearly 
interpolating the fitted values determined at the nearest simulated mass points. The fraction of events in 
the Voigtian part of the signal, /y, is ~ 60%(~ 70%) at 120 GeV, decreasing to ~ 58%(~ 55%) at 
400 GeV for the type-III seesaw signal (VLL). The uncertainties on the fit parameters of the signal p.d.f. 
are incorporated as Gaussian-distributed nuisance parameters in the combined fits to data. 

The combined background model consists of five different p.d.f.s, corresponding to ITZ, ZZ, Z + y, redu¬ 
cible, and the sum of the tt + V and triboson backgrounds. The leading ITZ and ZZ backgrounds are both 
modelled with a modified Bukin function [52], a three-parameter function designed to model peaks with 
asymmetric tails using the convolution of a Gaussian and an exponential function. To mitigate the impact 
of MC statistical uncertainties, the parameterizations for the Ai and SFonly categories are determined 
from the combination of all three categories; for the 37’-i-jj category, a separate parameterization is used 
to account for possible kinematic effects from the two additional jets. The uncertainty on the shape of 
the Am distribution predicted by the generator is taken into account with a Gaussian-distributed nuisance 
parameter multiplying a template given by the difference befween fhe p.d.f.s defermined from Sherpa and 
powHEG-Box. Finally, in fhe 3^-i-jj and 3Fonly cafegories, fhe rafio of normalizations of fhe WZ and ZZ 
backgrounds is fixed fo fhe prediction from MC simulation, due fo fhe inabilify of fhe fif fo resolve fhe 
similar shapes of fhe Am disfribufions. In fhe A£ cafegory, fhe confribufion from WZ evenfs is negligible, 
so only fhe ZZ background normalizafion is leff as a free paramefer. 

The mosf imporfanf remaining backgrounds are due fo reducible processes and Z + y production. The 
reducible background is parameferized wifh a Landau disfribufion, defermined from fhe dafa-driven es- 
fimafe described above. Due fo fhe low number of evenfs in fhe reducible background esfimafes, fhe 
corresponding shape is obfained by fiffing fhe combined disfribufion of all six signal regions. 

The Z{€{) + y background confribufes significanfly only fo fhe Z + e, 3Gonly cafegory. This background 
is modelled wifh fhe sum of a Landau disfribufion and a Gaussian funcfion. The normalizafions of fhe 
reducible and Z -i- y backgrounds are consfrained fo fhe expected values, each wifh a Gaussian-disfribufed 
uncertainly of 30%. 

Finally, fhe background confribufions from friboson and tt + WjZ producfion are modelled fogefher using 
a Landau disfribufion, due fo fheir similar shape and small confribufion fo fhe folal background eslimale. 
The normalizafion in fhe combined fif is also consfrained fo fhe expecfed value from simulation, wifh a 
30% Gaussian-disfribufed uncerlainly. 


8 Results 

The folal number of evenfs observed in each signal region is shown in fable 3, along wifh fhe eslimaled 
backgrounds before and affer fiffing fhe folal background model fo fhe dala. The corresponding Am dis- 
Iribufions for fhe pre-fif background estimates and fhe dafa are shown in figure 3. The signals expecfed 
for fhe VLL model wifh mL± - 140 GeV and fhe lype-III seesaw model wifh mL± = 300 GeV are super¬ 
imposed on fhe background as illuslralive examples. The dala agree wifh fhe background expeclalion in 
all cases, and no clear peak indicating resonanl Irileplon production is seen in any of fhe signal regions. 
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Figure 3: The Am - rnje - mt+f- distributions for the 4{ (top), Sf+jj (middle), and 3^-only (bottom) categories, 
divided into the Z + e (left) and Z + jj (right) flavour channels. The observed data are shown as black points, while 
the pre-fit background expectations are shown in the coloured histograms. Also shown are examples for signal 
contributions for a 140 GeV L* in the VLL model and a 300 GeV L* in the type-III seesaw model. The error bars 
on the data points represent statistical uncertainties, and the shaded bands represent the systematic uncertainties on 
the background predictions. 
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Table 3: Observed and expected number of events in the six signal regions, before and after the combined unbinned 
maximum-likelihood fit. The pre-fit uncertainties represent the total systematic uncertainties on the background 
estimates. The post-fit uncertainties are determined by the maximum-likelihood fit. 



'Zj + 6 

Z ± p 

Process 

41 SR 

3/±77 SR 

3/-only SR 

4/ SR 

3Z± 

././SR 

3Z-only SR 


Before combined background-only fit 

ZZ 

10.9 ±0.6 

11.7 

±0.8 

91 

±5 

21.4 ±1.1 

7.5 

±0.6 

90 

±5 

WZ 

0.08 ±0.01 

35.3 

±3.1 

337 

±28 

— 

46 

±4 

480 

±40 

Z -1- 7 

— 

2.3 

±0.8 

35 

± 11 

— 





Reducible 

— 

1.6 

±0.5 

38 

± 14 

1.5 ±0.3 

8.8 

±3.0 

79 

±22 

ft + V, VV 

1.2 ±0.2 

7.8 

± 1.7 

2.3 

±0.4 

1.5 ±0.2 

9.5 

±2.1 

3.3 

±0.5 

Total Background 

12.2 ±0.7 

59 

±4 

504 

±34 

24.4 ± 1.2 

72 

±6 

650 

±50 


After combined background-only fit 

ZZ 

15 ±4 

13.4 

±2.3 

107 

±9 

22 ±5 

10.1 

± 1.6 

88 

±8 

WZ 

0.08 ± 0.03 

39 

±6 

393 

±28 

0.02 ± 0.02 

56 

±9 

460 

±40 

Z + y 

— 

2.2 

±0.8 

34 

± 11 

— 





Reducible 

— 

1.8 

± 1.2 

37 

± 13 

1.8 ±0.9 

10.2 

±2.8 

92 

±24 

ft + V, 

1.1 ±0.3 

7.5 

± 1.7 

2.5 

±0.6 

1.5 ±0.4 

9.1 

±2.1 

3.3 

±0.8 

Total Background 

16 ±4 

64 

±7 

574 

±34 

25 ±5 

85 

± 10 

640 

±40 

Data 

16 

64 

573 

25 

86 

651 


Good agreement is seen between the pre-fit and post-fit normalizations for the Ai and 3^-i-jj categories 
in the Z + n flavour channel. The largest change in normalization due to the fit is in the A£ category for 
the Z + e flavour channel, where the fitted ZZ normalization exceeds the prediction by 35%. The WZ 
and ZZ normalizations increase by roughly 15% in the 3{+]] and 3^-only categories in the Z + e flavour 
channel, and 30% in the 3^-i-jj category in the Z + ji flavour channel. The projections of the fit results in 
the background-only hypothesis are shown in figure 4 for the combination of the three categories in each 
flavour channel. 

The data are well described by the combined fit to the three categories in each flavour channel. The con¬ 
sistency of the data with the background-only hypothesis is evaluated by scanning the local po'Value for 
the Am distribution in 3 GeV intervals for signal mass hypotheses in the range 100-400 GeV for the VLL 
model, and 100 - 500 GeV for the seesaw model, using the unbinned maximum-likelihood fit described 
in section 7 with the signal strength set to zero. The po'Value, which corresponds to the probability to 
observe at least as many events as observed in the present measurement assuming the background-only 
hypothesis, is calculated using the frequentist hypothesis test based on the profile likelihood ratio test 
statistic and approximated with asymptotic formulae [53]. The minimum po'Value is po = 0-02 at a mass 
of 183 GeV for the Z + e flavour channel, and po - 0-05 at a mass of 109 GeV for the Z -i- p flavour 
channel. 

Since no significant excess above the background expectation is observed, the fit model is used to derive 
95% confidence level (CL) exclusion limits on the heavy lepton pair-production cross section, cr, using the 
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Figure 4: Projections onto the Am variable of the background-only unbinned maximum-likelihood fits, shown su¬ 
perimposed on the data with the three categories in each flavour channel added together. The Z + e flavour channel 
is shown in (a), and the Z + /u channel is shown in (b). The contributions of the separate background components 
to the total background-only fit are also shown. The error bars on the data points represent statistical uncertainties. 
Good agreement is observed between the background model and the data. 


CLs method [54] . The limits are shown for the VLL model in figure 5, and for the type-III seesaw model in 
figure 6, evaluafed in fhe same 3 GeV infervals as fhe po'Values. The VLL model is excluded for elecfron- 
only mixing in fhe heavy lepfon mass ranges 129-144 GeV and 163-176 GeV, wifh an expected exclusion 
in the range 109-152 GeV. The corresponding observed (expected) exclusion for the muon-only mixing 
scenario is 114-153 GeV and 160-168 GeV (105-167 GeV). The significantly higher production cross 
sections for the type-III seesaw model lead to an observed (expected) exclusion in the electron-only 
mixing scenario in the heavy lepton mass range 100^30 GeV (100-436 GeV). For the muon-only mixing 
scenario, the observed exclusion is in the ranges 100-401 GeV and 419-468 GeV, while the expected 
exclusion is 100^19 GeV. 

The constraints shown in figures 5 and 6 are relevanf fo fhe specific VLL and fype-III seesaw models 
considered, and are nol necessarily applicable fo ofher scenarios predicfing frilepfon resonances wifh an 
infermediafe Z boson. A more model-independenf observable is fhe visible cross section, cTyis, defined as 
fhe number of observed evenfs with Z-i-Ginduced trilepton resonances for a given resonance mass divided 
by the integrated luminosity of the data sample, 20.3 fb“'. The 95% CL upper limits on cTyis, denoted 
are derived from a fit to each flavour channel with /y = 1, i.e. using only the peak component of 
the signal. The results for the two flavour channels, derived from the inclusive event selection without 
dividing the events into the three categories, are shown in figure 7. 

The limifs on cTyis can be used fo fesf specific models afler faking info accounf fhe model’s accepfance 
wifh respecf fo a fiducial volume, J?l, and reconsfrucfion and selecfion efficiency of evenfs wifhin fhe 
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Figure 5; 95% CL upper limits on the vector-like lepton cross section. The left (right) plot shows the limits assuming 
100% branching fraction to e/vg (ju/v^). The solid line shows the observed limit. The dashed line shows the median 
expected limit for a background-only hypothesis, with green and yellow bands indicating the expected fluctuations 
at the +lcr and +2cr levels. The limit is evaluated in 3 GeV intervals. 




Figure 6; 95% CL upper limits on the type-III seesaw production cross section. The left (right) plot shows the limits 
assuming 100% branching fraction to e/Ve (plVf^). The solid line shows the observed limit. The dashed line shows 
the median expected limit for a background-only hypothesis dataset, with green and yellow bands indicating the 
expected fluctuations at the +lcr and +2cr levels. The limit is evaluated in 3 GeV intervals. 
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fiducial volume, egd. The 95% CL upper limit on the cross section for the model is given by: 


(^95 


cr 


vis 

95 


J?1 X efid 


( 2 ) 


The acceptance ^ is defined as the probability for generated signal events to lie within a fiducial volume 
defined by the kinematics of the generated leptons. The leptons are considered at particle level, i.e. after 
parton shower and hadronization and with lifetimes longer than 10“^^ s, and are dressed, including the 
contributions from radiated photons within a cone of AR = 0.1. The fiducial volume requires that events 
contain an L- decaying to a prompt electron or muon and a Z boson that then decays to electrons or muons. 
The three leptons from the L- decay are required to have px > 15 GeV and lie within \r]\ < 2.5, with at 
least one lepton satisfying px > 26 GeV. Two of the leptons must form a same-flavour opposite-sign pair 
with a mass within 10 GeV of mz, and the Z boson and the off-Z lepton must be separated by AR < 3. 
The events are divided into flavour channels according to the flavour of the off-Z lepton. For the VLL 
and type-III seesaw models used in this analysis, the acceptance of events containing an Lr —> Z{(£)£ 
decay to fall within the fiducial volume is in the range 60%-65% for most of the mass range, decreasing 
at higher masses due to the cut on the AR between the Z boson and the off-Z lepton. The acceptance 
decreases at low masses due to the lepton px requirement, reaching 30%-35% at mx± = 100 GeV. 

For type-III seesaw and VLL events within the fiducial volume, egd ranges from 20% to 49% if the other 
heavy lepton decays to a neutrino and a W, Z, or H boson If the other heavy lepton decays to an electron 
or a muon, the efficiency is 10%-20% lower, due to the increased probability of incorrectly selecting the 
off-Z lepton. The event selection efficiencies for the type-III seesaw model in scenarios where the second 
heavy lepton decays to a fT boson are shown in figure 8 as a function of ; the efficiencies for scenarios 
where the second heavy lepton decays to a Z or // boson and for the VLL model are consistent with these 
efficiencies within the statistical uncertainties. 



Figure 7: Upper limits at 95% CL on cxvis for the Z -i- e (left) and Z + p (right) flavour channels, derived without 
dividing events into the three categories. The limits are evaluated in 3 GeV intervals. 


^ Note that the quoted efficiencies are dependent on the modelling of the polarization of the Z bosons, due to the requirements 
imposed on lepton isolation and separation. 
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Figure 8; Efficiencies for reconstructing and correctly identifying the V" decay in events within the 

fiducial volume for the type-III seesaw model. The left (right) plot shows the efficiencies for events containing a 
L* ^ Z{t€)e (L* ^ Z({£)iJ.) decay. The decay of the second heavy lepton is specified in the legend. The shaded 
bands show the statistical uncertainty. 


9 Conclusion 

A search for trilepton resonances decaying to a Z boson and an electron or a muon has been presented. 
The search is based on pp collision data taken at = 8 TeV corresponding to an integrated luminosity of 
20.3 fb“^, collected by the ATLAS experiment at the CERN Large Hadron Collider. Events are selected 
requiring at least three electrons or muons with high transverse momentum, with two of the leptons 
consistent with the decay of a Z boson. The events are categorized based on the presence or absence of 
additional leptons or dijet pairs in the event consistent with the decay products of a second heavy lepton, 
and separated into channels based on the flavour of the lepton associated with the Z boson to form a heavy 
lepton decay candidate. Using the difference between the trilepton and the Z boson candidate masses, a 
search for a narrow resonance is performed in each of these categories using a maximum-likelihood 
fit of parameterized signal and background shapes to the data. No significant excess above Standard 
Model predictions is observed, and 95% CL limits on the production of trilepton resonances beyond 
the Standard Model are derived. The results are interpreted in the context of two models of new heavy 
leptons decaying to three charged leptons. In the vector-like lepton model, new heavy charged leptons are 
excluded in the mass range 129-176 GeV (114-168 GeV) for electron-only (muon-only) mixing, except 
for the interval 144-163 GeV (153-160 GeV). In the type-III seesaw model, the corresponding exclusion 
is in the mass range 100-430 GeV (100^68 GeV) for electron-only (muon-only) mixing, except for the 
interval 401-419 GeV in the muon case. Limits are also set on the visible cross section of trilepton 
resonance productions, and fiducial efficiencies are derived to facilitate model testing. 
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